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Abstract. The propagation of UHECR nuclei for A = 1 (protons) to A = 56 (iron) from cosmo- 
^"J ' logical sources through extragalactic space is discussed in the first lecture. This is followed in the 

f"| , second and third lectures by a consideration of the generation and propagation of secondary parti- 

Q-i cles produced via the UHECR loss interactions. In the second lecture we focus on the generation 

q \ of the diffuse cosmogenic UHE-neutrino flux. In the third lecture we investigate the arriving flux 

j_i . of UHE-photon flux at Earth. In the final lecture the results of the previous lectures are put to- 

' gether in order to provide new insights into UHECR sources. The first of these providing a means 

C$ [ with which to investigate the local population of UHECR sources through the measurement of the 

1—1 i UHECR spectrum and their photon fraction at Earth. The second of these providing contraints on 

the UHECR source radiation fields through the possible observation at Earth of UHECR nuclei. 
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LECTURE 1: UHECR PROTON AND NUCLEI PROPAGATION 

(N ■ 
O 

Q\ | Over the last century, Cosmic Rays (CR) have been detected arriving at Earth. The 

energy range with which these particles arrive spans over ten decades in energy from 
10 9 eV to 10 20 eV. For a large part of this energy range, the CR spectrum is found 
to fit a polynomial description of the form dN/dE E~ a , with a the spectral index, 
taking value of roughly 2.7 up to 10 15 eV (1 particle m 2 yr _1 ), changing above this 
energy to a = 3.1 (producing a feature called the "knee" [1]). Further on, at higher 
energies, a second "knee" in the spectrum is observed near 10 17 eV. Following this, 
at energies just above 10 18 eV (1 particle km 2 yr _1 ), another sudden change occurs 
in the spectrum, with it apparently returning to a spectral index value of a = 2.7 (the 
"ankle"). CR with energies over 10 18 eV are referred to as Ultra High Energy Cosmic 
Rays (UHECR). From a simple consideration of the Larmor radius of CR in the Galaxy it 
is seen that UHECR containment within the Galaxy's /lG magnetic fields is problematic, 
most probably requiring their sources to be extragalactic in origin at these energies. 
Being charged particles, UHECR will execute trajectories with non-negligible curvature 
in nG extragalactic magnetic fields present on distance scales >Mpc from their sources 
[2]. This would lead to a washing out of their directional (source) information on such 



scales. 

With direct CR detection only being possible in the low energy region < 10 14 eV, 
where the arriving fluxes are sufficiently large to allow satellite and balloon spec- 
troscopy, a direct handle on the CR composition is only possible for these lower energy 
particles. At higher energies only indirect measurements are viable through the analysis 
of the profile and content of the particle shower created in the atmosphere. However, for 
the hadronic part of the shower, the composition of the primary particle plays an impor- 
tant role in the shower development, dictating the eventual shower profile and content. 
Thus, the composition information is encoded in the CR showers detected. 

In this first lecture we consider the propagation of UHECR nuclei for A = 1 (protons) 
to A = 56 (iron) from cosmological sources through interstellar space to Earth. Firstly, 
the physics governing the propagation of UHE-protons will be discussed. This is fol- 
lowed in the second section by a discussion of nuclei UHECR propagation. In the third 
section we study the UHECR spectra observed at Earth for different primary particles 
injected at the sources. Finally, we end the lecture with a brief description of an analytic 
approach for UHECR nuclei propagation. 

1.2 COMPOSITION OF COSMIC RAYS 

As already mentioned, the flux of UHECR is sufficiently low that a present-day satellite 
born flux measurement, if attempted, would need to wait longer than ~ 10 yrs just 
to detect one. Instead, in order to detect a significant amount of events it is necessary 
to use the Earth's atmosphere as the detector, as employed by current day UHECR air 
shower detectors such as the Pierre Auger Observatory [3] (Auger) in Argentina. In 
these experiments, indirect detection of UHECR is achieved via the secondary particles 
generated in a particle shower created through the interaction of a (primary) CR with the 
nuclei of atoms in the atmosphere. The price payed through this approach, with present 
uncertainty in the description of the hadronic shower development, is the loss of precise 
composition information of the primary particle. Nevertheless, there are some particle 
shower characteristics which do intimately relate to the primary particle composition, 
they are: 

• the distance from the "top of the atmosphere", measured in column depth, at which 
the number of particles in the shower is maximum (X max ), X^ ax (E) > X^l x (E) 

• the number of muons in the shower (N^), N^ e (E) > Nj^(E) 

We here focus on the first of these composition indicators. The energy of the shower 
particles energy falls below a critical value E CT [ t . Through a Heitler model [4] descrip- 
tion for the electromagnetic shower, with primary energy Eq, the conservation of energy 
in the shower requires, Eq = E CT n 2" crit ■, where n CI [ t is the number of particles in the 
shower when the particles in the shower have energy ^crit.- This X max = n cr i t .A ln(2) = 
X ln(Eo/E CI i t ) where A is the mean free path. Typically E CT [ t = 85 MeV for electromag- 
natic showers (f/e). 

Fig. 3 in ref. [5] demonstates that recent Auger X max measurements presently support 
arguments for a mixed composition of UHECR arriving at Earth. With such measure- 
ments as motivation, along with UHECR proton propagation we also review the main 



aspects of UHECR nuclei propagation in the following sections. Indeed, it is worth stat- 
ing that the CR composition observed at Earth may be quite different from that at injec- 
tion. Intermediate mass or heavy nuclei injected in a distant CR acceleration will gradu- 
ally disintegrate into lighter nuclei and nucleons as they propagate through intergalactic 
space [2]. 



1.3 UHECR PROTON ENERGY LOSS INTERACTIONS 

The photon targets for pion production are given by the Cosmic microwave background 
(CMB) and the Cosmic infrared background (CIR). The specific number density of the 
CMB photon spectrum constitutes the main target for py interactions and is described 
by the Planck's law for a black body of temperature kT = 2.3 x 10~ 4 eV (2.7 K). 
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where c is the velocity of light in a vacuum and h is Planck's constant. 

Over large enough distance scales through this background, UHECR protons will 
undergo interactions with these photons leading to the production of electron positron 
pairs (p + Y^ p + e + e~) and pions (p + Y^ p/n + n° 

The attenuation rate calculation for both these (py) proton energy loss processes may 
be expressed as 
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where Y = E p /m p , m p is the protons mass, e r is the photon energy, o pr is the pY 
interaction cross section and K p is the inelasticity of the proton for such interactions 
(Kp = AE p /E p ). 

Pair production- In the rest frame of the proton, the photon threshold energy for pair 
1 MeV. The inelasticities of these interactions typically going as 
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Arn^c 4 /(m p c eL), where eL is the colliding photon's energy in the proton's rest 
frame and m e is the electron rest mass [6]. 

Pion production-In the rest frame of the proton, the photon threshold energy for pion 
production is e' th ~ 145 MeV. The inelasticities of these interactions, for low pion 

multiplicities, go approximately as K p ~ (m^. + 2m p e'y)/2(mp + 2m p eL), where e' r is 
the colliding photon's energy in the proton rest frame, and m n is the pion rest mass [7]. 

For an approximate description of R(E p ) for pion losses, we may describe the cross- 
section as a top-hat function, 
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where Oa ~ 0.5 mb is the peak value of the cross section (A resonance), which occurs at 
e' A « 340 MeV, and 8 « 100 MeV [8]. With inelasticities of pion production processes 
occurring close to threshold typically taking values, K p ps 0.2, 
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where Iq = 5 Mpc and x = 10 2053 eV/E p (note- this calculation assumes only the 
presence of the CMB photon target exists). 

The rate curves for these processes, shown in the left-panel of fig. 1, may now be well 
understood, with their maximum values being of the order, 



E —— for pair production 

for pion production. 



1.4 COSMIC RAY PROTON PROPAGATION 



Once UHECRs have been accelerated and escaped their source region, their journey 
through extragalactic space begins. UHECR protons, of energies above 10 20 eV, are not 
expected to be deflected sufficiently by Galactic or extragalactic magnetic fields that 
their source directional information will be removed. Along with this, their propagation 
is limited through their inelastic collisions with cosmic background photons, attenuating 
their propagation to < 30 Mpc distance scales [9, 10] (see eqn 4), i.e. about the size of 
the local supercluster of galaxies. Recent measurements by the Auger detector indicate 
that both that the arriving UHECR spectrum is not isotropic at high energies and that 
a suppression is observed in the arriving UHECR flux [11, 12]. However, more data is 
required for further clarity as to the cause of these results. In preparation for developing 
such an understanding, we here outline the theoretical tools required in order that 
calculations of UHECR proton propagation may be carried out. 

With very small inelasticities for the proton in each interaction, to a very good 
approximation pair production can be treated as a continuous energy loss process. Photo- 
pion production interactions, however, involve very large inelasticities, so may not be 
treated as a continuous loss process and it is necessary to use Monte Carlo techniques 
[13]. 

It is also possible for UHECR protons to produce pions through interactions with 
the cosmic infrared background [2]. Although the rates for these interactions are sub- 
dominant in comparison to energy losses from pair production, they can be important 
in determining the spectrum of cosmogenic UHE-neutrinos produced through UHECR 
proton propagation [14]. 

At present there are conflicting claims concerning the existence of the GZK cut-off, 
with two of the three experiments with largest total exposure reporting a cut-off signature 
[15, 16, 17]. Interestingly, departures from the "vanilla-sky" expectation for a GZK 
suppression are not so difficult to obtain with both exotic new physics and less-exotic 
scenerios through alterations to some of the standard assumptions. Indeed, with regard 
these less exotic scenarios, the actual feature to be expected for the cut-off is dependent 



on several factors such as the local source distribution, composition, cross-over energy 
at which the extragalactic UHECR flux dominates [18]. In particular, if UHECRs are not 
protons but consist of a substantial amount of heavy nuclei (see for example [19, 20]) 
then the spectrum will be altered from the usual expectation. Although evidence for the 
presence of heavy nuclei in UHECRs has been around for some time, the information 
available is still imprecise. 



1.5 AN ANALYTIC DESCRIPTION FOR PROTON PROPAGATION 

For a source at distance L, the ratio of the number of protons, N n (E p ,L), arriving with 
energy E p having undergone n pion production interactions to the initial population, 
Nq(E,0), of protons of energy E is [18] 

^^ = £ ;oC -.exp(-L/ W n-l r , (6) 

^vo^uj m=Q p =0 l m l p 

with 
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where K p is the inelasticity, Iq = 1 Mpc and x = 10 205 eY/E. 

Thus, the distribution of protons of a given energy having not undergone any pion pro- 
duction interactions is given by a simple exponential decay function. Similarly, the dis- 
tribution of mono-energetic protons having undergone one pion production interaction 
but not having undergone a second pion production interaction is given by the difference 
of two exponential functions. 

A demonstration of the success of this description may be found in [18], for which 
the arriving UHECR flux from the local region is obtained and compared with the 
numerical result. Interestingly, a similar expression as to that provided above (eqn 6) 
may also describe the distribution of nuclear species produced following UHECR nuclei 
propagation, this will be discussed further in section 1.7. 



1.6 COSMIC RAY NUCLEI PROPAGATION 

In this section we study the intergalactic propagation of UHECR nuclei and their spec- 
trum at Earth for different species injected at source. Nuclei undergo different interac- 
tions with the CMB/CIB 

( A z)y^ H A ^ e+e ~ forl ° 19 - 7 eV < E (A,Z) < 1020 2 eV m 
1 ' )Y \{A- l,Z)n/p for£ (A Z) < 10 19 7 eV and E {A)Z) > 10 202 eV 

Moreover, nuclei suffer photon-disintegration 
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Figure 1. Left-Panel: Energy loss rates due to proton interactions with the CMB background photons. 
Right-Panel: Energy loss rates due to nuclei interactions with CMB and CIB background photons. 
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where T(a.z) = E(A,z)/Am p c 2 , A and Z denote the total number of nucleons (pro- 
tons+neutrons) and charge of the nucleus, i„ and i n are the number of protons and neu- 



tron broken off from a nucleus interaction and OA.z,i p ,i„y(£y) is the cross-section (see [2]). 
In the right-panel of Fig. 1 we present the energy loss lengths for iron nuclei through pair 
creation photo-disintegration interactions. 

In general, during propagation, a heavy nuclei will undergo many photo-disintegration 
reactions, cascading down in atomic number and charge, generating secondary UHECR 
protons, neutrons and alpha particles along the way, each of which may also contribute 
to the UHECR spectrum at Earth. The amount of photo-disintegration depend on the 
amount of time they spend in the radiation fields, the strengths of the radiation fields 
that they propagate through, and the composition of the CRs leaving the accelerating 
region (source) [2] . 

For the purpose of CR propagation simulations, for energies above 10 19 eV, only 
sources up to redshift of 0.5 really contribute to the CR flux at Earth, this limit being set 
by pair production energy losses of CRs during propagation. In general terms, during 
the past decade there has been quite a bit of work on the propagation of UHECR nuclei 
(see for example ref.s [21, 2]). Usually, the complex process of the photo-disintegration 
of UHECR nuclei into lighter nuclei and nucleons has been so far addressed using 
Monte Carlo techniques. More recently, a new analytic approach has been implemented 
with very good results [22]. We will describe this subject briefly as a final point in our 
discussion. 



1.7 AN ANALYTIC DESCRIPTION FOR NUCLEI PROPAGATION 



Denoting the number of nuclei with atomic number A by Na, the differential equation 
describing the population of a nuclei state is given by 

dN A + Na + Na + = N A +i + N A +2 + 9 
dL Ia^a-i Ia+\^a Ia+2-^a 

where L is the distance traveled and l^j is the interaction length for state i to disintegrate 
in state j. In order to simplify the calculations, we consider the case in which only 
a single nucleon loss process occur. This reduces the number of states in the system 
dramatically, along with the number of possible transitions. The equation (9) simplifies 
to 

dNA + NA = NA+l 

dL I a Ia+i 

The solution to this set of coupled differential equations, constrained by the initial 
conditions N n (L = 0) ^ and N n (L = 0) = for A n is given by 

^ = £ W --exp(-f) n t^t ai) 

^n\y) m=A hn p=A (^ m ) Lm P 

A demonstration of the success of this description may be found in [22], for which 
the arriving UHECR nuclei flux and composition is obtained and compared with the 
numerical result. 



1.8 CONCLUSION 

In this lecture, the dominant energy loss interactions of UHECR protons during their 
propagation through extragalactic radiation fields were described. An analytic descrip- 
tion for proton propagation through such radiation fields was also provided. Following 
a motivation for the consideration of UHECR nuclei propagation, from measurements 
made by present generation UHECR detectors such as Auger, we also considered the 
dominant energy loss in interactions of UHECR nuclei. An analytic description for this 
situation was also provided. 



LECTURE 2: THE COSMOGENIC UHE-NEUTRINO FLUX 

During the propagation of ultra-high energy (UHE) protons through extragalactic space, 
a finite but small probability exists that they may undergo an interaction with each of 
the background photons they encounter. The most abundant of the background photons 
encountered, with sufficient energy to partake in these interactions, existing in the 2.7 K 
microwave (10~ 3 eV) and infrared (10~ 2 eV) backgrounds. The large distances present 
in extragalactic space, however, provides sufficient column densities of such photons 
between the source and Earth to make such interactions likely, leading to the limitation 



of UHECR proton propagation in extragalactic space (assuming rectilinear propagation) 
to distances less than t\l e ~x\ , where x = 10 2053 eV/E p , E p is the proton's energy, 
and Iq is 5 Mpc. The energy fluxes of secondaries produced through such interactions 
(p + 7 ~^ p / n + n° / n + ) leads to the generation of both neutrino and photon fluxes. 
However, we here will focus on the production of the (cosmogenic) UHE-neutrinos, 
leaving the photon flux for the following lecture. 

In the next section we will focus on the photo-pion producing interactions. This 
will be followed by a section describing the expected (cosmogenic) UHE-neutrino 
flux produced by UHECR propagation through the cosmic background radiation fields, 
indicating how it is dependent on the CR composition at the source. 

2.1 SOURCE DISTRIBUTIONS 

In order to calculate the cosmogenic UHE-neutrino spectrum arriving at Earth it is 
necessary to assume an UHECR energy spectrum injected at source, dN p /dE p \ & r , 
and the distribution of sources, described by the number of sources in shells of redshift 
bins redshift surrounding the Earth, dN sowce /dz. 

The energy spectrum assumed to be produced by our sources are motivated by Fermi 
first order non relativistic shock acceleration theory. In this theory, particles gain energy 
through consecutive crossings of a non-relativistic shock wave. The energy gain per 
crossing is found to be AE = ^[3Eq, where /3 is the relative velocity of the two plasma's 
either side of the shock in units of c and Eq is the initial energy of the particle. At 
the same time as the energy of the particles increases with each crossing, the number 
of particles taking part is expected to decrease, with particles being lost from the 
system through advection (being washed downstream). The number of particles lost 
in each crossing being given by, AN = — |j8iVb, where N Pj o is the initial number of 
particles in the system. After n crossings, the energy of the particles will be E = 
(1 + |/3 )"£(), whereas the number of particles will be N = (1 — |/3)'Wo. From these 
distributions, one sees that of order n ~ l/j8 encounters are needed to significantly 
modify the particle energy and the number of particles. For /3 1, the altered particle 
distribution becomes dN/dE °c E~ 2 , the spectrum of Fermi accelerated particles, with 
spectral index a = 2. Such a theoretical E~ 2 spectrum is consistent with the inferred 
electron/proton spectra from y-ray observations in the GeV-TeV domain, produced by 
CR electron/proton energy loss processes close to the source [23, 24, 25]. 

From this motivation, we assume a UHECR spectrum at the source described by a 
power law with an exponential cutoff dN/dE\ source oc £-« g -£/£max i where the maxi- 
mum energy is 1 E max = (Z/26) x 10 22 eV. 

For the distribution of sources per co-moving volume the quasar luminosity density 
evolution[26] is adopted 



following Hillas criterion type arguments, the maximum energy attainable by a source is expected to be 
proportional to the charge Z of the particle being accelerated 
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2.2. COSMOGENIC UHE-NEUTRINO PRODUCTION 



Neutrinos created during UHECR propagation are known as cosmogenic UHE- 
neutrinos. As discussed in the first lecture, after leaving the source UHECR protons 
may interact with the CMB photons creating charged pions through photo-pion pro- 
duction interactions (see section 1.3 in the first lecture). Following their decay, charged 
pions produce a flux of neutrinos, 

tu + ^v^+fi + ^ v M +e + + V e + v M . (13) 
Neutrinos are also produced via beta decay of secondary neutrons 

n^p + e~ + V e . (14) 

The energy of the neutrinos resulting from py reactions is E v ~ 0.05£p, and in the case 
of neutron decay, E v ~ 0.0005£p. 

Following the assumption about the energy and spatial distributions of the sources, 
discussed in the previous section, we here obtain the subsequent cosmogenic UHE- 
neutrino flux produced. In Fig. 2, the charged pion decay generated neutrino flux at 
its maximum (the higher energy peak) is roughly three times the value of the neutron 
decay generated electron neutrino flux (lower energy peak). The normalisation of the 
cosmogenic UHE-neutrino energy flux, with a value ~ 10 eV cm' 2 s _1 sr _1 , is compa- 
rable to the UHECR energy flux at 10 19 eV, as might be expected by the large decrease 
in the attenuation length of UHECR protons at slightly higher energies, whose energy 
flux feeds into the UHE-neutrino flux. 

Most models for the production of cosmogenic UHE-neutrinos assume that UHECRs 
are protons. This leads to an observable flux of neutrinos in detectors such as IceCube 
or KM3net. However, the existence of arriving UHECR nuclei can imply that these 
cosmogenic UHE-neutrino fluxes are suppressed by factors as large as ~ 100 (see [27]), 
making it hard to detect the (previously considered) "guaranteed" flux of cosmogenic 
UHE-neutrinos. Instead, we here calculate also the neutrino flux expected for the case 
of UHECR nuclei. 

During propagation, UHECR nuclei experience photo-disintegration into their con- 
stituent nucleons in scattering off the CMB and CIB photons [2]. The interaction of 
the daughter nucleons with these background photons may create charged pions which 
will go on to decay producing a cosmogenic UHE-flux of neutrinos. In the nuclei rest 
frame, the photo-disintegration cross section peaks at a photon energy of ~ 30 MeV (gi- 
ant dipole resonance). Since in the lab frame, the average CMB-CIB photon energy is 
~ 10 eV, the nucleus must have an energy of approximately 10 20 eV for the interac- 
tion to take place. This energy is roughly the same as that needed for protons to interact 
with the CMB photons. Thus, the GZK suppression for nuclei can be expected to occur 
at similar energies to that for protons. 




Figure 2. Left-Panel: The calculated cosmogenic UHE-neutrino flux produced by the cosmic popula- 
tion of UHECR protons assuming the energy spectrum and spatial distribution given in section 2.1. This 
result sits in good agreement with that calculated in ref. [28]. Right-Panel: The calculated cosmogenic 
UHE-neutrino flux produced by a cosmic population of UHECR protons and nuclei assuming the energy 
spectrum and spatial distribution of sources given in section 2.1. This result was obtained in ref. [27]. 



For a given spectrum and composition at the sources, the calculated spectrum at 
Earth must be consistent both with the observed UHECR spectrum and the X max (the 
atmospheric depth of the shower maximum) data. Whereas, the X max data seems to 
find better agreement for an all-iron spectrum compared to an all-proton composition 
at source (see for instance [27]), the UHECR spectrum does not give a clear indication. 

The cosmogenic UHE-neutrino flux for CR nuclei is presented in the right-panel 
of Fig. 2, together with the all-proton "guaranteed" neutrino flux (dashed line). The 
upper and lower solid line corresponds, respectively, to the highest (dominant proton 
composition) and lowest (dominant iron admixture) neutrino flux compatible with the 
data. It can be seen that the lowest possible value is two orders of magnitude below 
the previously considered all-proton cosmogenic UHE-neutrino flux. If this is the case, 
it will not be possible to observe the cosmogenic UHE-neutrinos in detectors such as 
IceCube or KM3net. 

2.3 CONCLUSIONS 

In this lecture, a motivation for injection spectrum spectral index close to 2 was provided 
through the consideration of the Fermi first order non-relativistic shock acceleration 
mechanism. Through simulations of sources with such spectra, with an assumed cos- 
mological evolution history, the diffuse cosmogenic UHE-neutrino flux was obtained. 
Further to this, with the Auger X max data presently suggesting that UHECRs inject a 
component of heavy nuclei at source, this diffuse cosmogenic UHE-neutrino flux cal- 
culation has also been carried out for injection spectra and composition consistent with 
both the arriving spectrum and X max Auger measurements. In this way, the cosmogenic 
UHE-neutrino flux calculation has been demonstrated to rest upon important underly- 
ing assumptions, with the flux typically determined being by no means guaranteed. The 
(cosmogenic) UHE-neutrino flux may be reduced up to a factor of 100 if instead of the 



all-proton (previously studied) scenario, a large heavy nuclei component at the sources 
exists [27]. 

LECTURE 3: THE UHECR PHOTON FRACTION 

3.1 UHE-PHOTON PRODUCTION THROUGH UHECR LOSSES 

Following the acceleration of UHECR by their sources, and subsequent escape, this flux 
of particles (if protons) may be attenuated by photo-pion production interactions (p + 
y ^ p/n + n°/n + ) with the background radiation photons. Through this attenuation, 
a subsequent secondary flux of UHE-photons is generated via the decay of neutral 
pions (7T° — > 77). From isospin considerations for the Delta-resonance decay, neutral 
pions are twice as likely to be produced as charged pions. However, close to threshold 
{£y ~ 145 MeV) where direct pion production dominates, charge pion production is 
more likely. Nevertheless, unlike high energy neutrinos produced through the decay of 
charged pions (7T+ /iz~ — > + v^/v^ — > e + /e~ + v M + v M + v e /v e ), these UHE- 

photons are not free to propagate through the Universe over cosmological distances. The 
UHE-photon flux is in fact attenuated through interactions with the background radiation 
field via pair production 77^ e + e~ on Mpc size scales, similar to the scales from their 
sources on which these photon fluxes are produced. 

UHE-electrons produced through such pair creation interactions may either subse- 
quently undergo inverse Compton cooling interactions off the background radiation field 
through high center-of-mass ey collisions, leading to a repeating cycle of pair produc- 
tion and inverse Compton scattering and the development of an electromagnetic cascade 
[29, 30, 31], or synchrotron cool in the extra-galactic magnetic fields [31, 32]. This is 
investigated in detail in the following section. To conclude this section we would like 
to remark that high energy photon production is an inevitable consequence of the GZK 
cut-off's existence. 

3.2 THE PROPAGATION OF UHE-PHOTONS 

In this section we discuss how the UHE-photon flux, produced near the UHE-proton 
sources, is attenuated by their interaction with different cosmic background radiations 
and extragalactic magnetic fields. UHE-photon interaction with background radiation 
fields may lead to an electromagnetic cascade generated before these generated UHE- 
photons reach Earth. Such a cascade being generated by the repetition of two processes: 
(1) pair creation through UHE-photon interactions with background photons and (2) 
photon production by UHE-electrons inverse Compton scattering background photons. 
However, with regards the second process, the UHE-electrons may also synchrotron cool 
in the extragalactic magnetic field. In the following, we shall look more closely into the 
physics of pair creation and investigate where the energy content of the electrons and 
positrons produced flows to next and how the dominant cooling channel is decided upon. 

In the center-of-mass frame the electron/positron pairs are produced with equal en- 
ergy. However, following a boost back to the lab frame, one of the electron's tends to 



take nearly all the energy. The average ratio between the interacting particle E e and the 
particle produced Ey, for s 3> 1, being approximately, 



TT ~ — -, f(s) (15) 

Ey As 

where s is the squared center-of-mass energy of the scattering process in units of 
(2m e c 2 ) 2 (for convenience since this sets s = 1 as threshold for the process), given by 

AE v E v g 

(16) 



(2m e c 2 ) 2 

for a head-on collision in the lab frame. A description of f(s) is shown in the right-panel 
of Fig. 3. 

For large s, one of the electrons produced carries away most of the original photon 
energy. What happens with this electron energy? The UHE-electrons may lose energy by 
the interaction either with the CMB radiation field (through Inverse Compton scattering) 
or with the extragalactic magnetic field (through the synchrotron channel). Since the 
energy density of both fields are Uy MB = 0.25 eV cm~ 3 and Ub = 10~ 8 eV cm -3 

(B = 3 x 10 -10 G), one might naively expect that the dominating interaction will be 
the field with the larger energy density (ie. the inverse Compton process). However, we 
here pay a closer look at the physics of these competing cooling channels. 

In the Thomson approximation, the center-of-mass frame is the electron's rest frame, 
so the background photon in the electron's frame simply has its momentum reversed 

during the interaction. This occurs when T e Ey g < m e c 2 , where T e is the electron's 
Lorentz factor and m e is its rest mass. For synchrotron photons, we consider the magnetic 
field as a virtual field of photons, with characteristic energies, 

E r g = (it-) m ec 2 , (IV) 



B 



cnt. 



where B cxlL = 4 x 10 13 G. Hence, with B ^nG a value for the virtual photons of Ey g « 
1(T 18 eV is obtained. For E e = 10 19 eV (T e = 2 x 10 13 ), the boosted virtual photon's 
energy in the electron's rest frame 10 eV (ie. well within the Thomson 

regime). 

On the other hand, for the CMB photons with Ey MB ~ 10~ 3 eV, such photons in the 

electron's rest frame (for the same energy electron) have an energy T e Ey g = 2 x 10 10 eV, 
thus the Thomson approximation is not applicable so the Klein-Nishina description (ap- 
plicable for T e Ey g > m e c 2 ) must be applied instead. Once in the Klein-Nishina regime, 
the center-of-mass frame can no longer be taken as the electron rest frame and the inverse 
Compton scattering cross-section is suppressed, leaving the dominant cooling channel 
less obvious. If the CMB is the only background radiation field and the extragalac- 
tic magnetic field is roughly nG in magnitude, the magnetic field (synchrotron) cool- 
ing channel wins and UHE-photons will undergo cooling interactions predominantly 
through the synchrotron channel. 
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Figure 3. Left-Panel: The interaction lengths of UHE-photons propagating through the background 
radiation fields in extragalactic space. Right-Panel: The function f(s) for yy interactions for an isotropic 
distribution of target photon momenta in the lab frame. 



However, if there exists a stronger cosmic radio background (CRB) component 
(Ey g « 1(T 6 eV) with £/Jr RB « 1(T 8 eV cm~ 3 then F e Ey g = 2 x 10 5 eV thus the Thom- 
son approximation still roughly holds and the inverse Compton channel can dominate. 
Under this condition, along with extragalactic magnetic fields of ~ 10~ 2 nG, the inverse 
Compton channel will dominate the electron/positron cooling. In the scenario with dom- 
inant inverse Compton cooling, an e/y cascade is produced. If this cascade developes 
with interactions occurring in the Klein-Nishina regime, the high energy particle simply 
changes from a neutral to charge state and back spending roughly 1/3 of it's time in the 
neutral (photon) state and 2/3 in the charged (electron/positron) state. 



3.3 THE PROTON/PHOTON ATTENUATION RATES AND RATIO 

AT EARTH 



In the last section we discussed the difficulties that a high energy photon undergoes 
in the Universe. We here look at a comparison of the proton and photon propagation 
through extragalactic space in order to study the photon/proton ratio expected at Earth. 
Considering UHE-photons as UHECR proton interaction products, it is necessary to 
calculate the UHE-photon interaction lengths with different possible realisations of the 
background radiation field field, as shown in the left-panel of Fig. 3. 

A clear difference in the interaction lengths of UHE-photons is seen to be expected for 
the Clark radio background value [33] and that provided purely by the radio component 
of the CMB, for UHE-photons with energies > 10 19 eV. Furthermore, a clear difference 
is seen to be expected between the interactions rates of UHE-photons in the Clark and 
Protheroe and Biermann [34] backgrounds for energies > 10 20 eV. However, for these 
calculations, the Clark radio background will be assumed. 

With the physics of UHE-photon interactions discussed in the previous section, the 
question as to what fraction of UHECR are photons may be addressed. We here discuss 
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Figure 4. Left-Panel: The UHE-photon flux, produced through UHECR K° interactions, arriving at 
Earth, the different colour lines in the figure depicting a breakdown in the contributions to the photon 
flux from different shell radii. Right-Panel: The photon fraction of UHECR arriving at Earth, the different 
colour lines in the figure depicting a breakdown in the contributions to the photon flux from different shell 
radii. 



the contribution from different source regions by separating out the fluxes produced 
from sources within the local — 240 Mpc region surrounding the Earth. We assume that 
UHECR sources have a homogeneous density distribution locally, with the number of 
sources scaling simply with the co-moving volume element under consideration. Along 
with this, we assume each source is produces UHECR with a power law energy spectrum 
and exponential cut-off. These source spatial and energy distributions being given by 

— oc 1 +Z 3 ^ E -a e -E p jE p (lg) 

dV c dE p 1 

with a = 2 and E^ 3 * = 10 20 5 eV. Following these assumptions, the arriving proton and 
photon fluxes were calculated for different sources shells of radii — 3 Mpc, 3 — 9 Mpc, 
9-27 Mpc, 27-81 Mpc, 81 -243 Mpc and energy range 10 19 > E p > 10 21 eV. The 
resulting arriving UHE-photon flux and photon fraction of UHECR are shown in the left 
and right panels of Fig 4. 

With searches for UHE-photons in UHECR data providing photon fraction upper lim- 
its [35], we here look into the photon fraction values we expect through UHECR energy 
losses. Investigations into the photon flux produced through the GZK interactions have 
demonstrated that significant photon fractions can be expected for particular assump- 
tions about the UHECR sources and extragalactic environment. In the "Clark radio, 
single" scenario (assuming the Clark radio background and extragalactic magnetic field 
stronger than 0.1 nG [33]) the arriving photon fraction is predominantly contributed 
from sources in relatively local shells, with distances tens of Mpc away. For the ho- 
mogeneous source distribution case, photon fractions between 10 3 and 10 2 can be 
expected with the photon fraction peaking at 10~ 2 , approximately the maximum energy 
the UHECRs are injected at (E™ x = 10 20 5 eV). The effect of changing this maximum 
energy are shown in Fig. 5 of the last lecture. Details about photon fraction can be found 
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Figure 5. Left: The UHECR proton flux arriving from a homogeneous distribution of sources for 
different values of E max . Right: The corresponding photon fraction of the arriving UHECR obtained from 
a homogeneous distribution of sources for different values of E m:iK . For this plot, the label "Clark radio, 
single" refers to the assumptions about the radio background and magnetic field strengths, for which 
we have assumed the Clark et al. [33] radio background and extra-galactic magnetic fields stronger than 
0. 1 nG. The injected proton spectrum assumed for these results had a spectral index of 0£=2. 



in [18, 32, 36] 



3.4 CONCLUSION 

In this lecture, the energy losses of UHE-photons in both magnetic fields and radiation 
fields have been discussed. We have demonstrated that this photon fraction provides 
a complementary measurement to that of the cut-off feature in the UHECR spectrum, 
allowing differentiation between the "tired source" and the GZK cut-off scenarios as 
well as carrying information about the energy injection spectrum at the source. Though 
measurements by Auger of the arriving UHECR flux photon fraction are only able to 
presently provide lower limit values (~ 10~ 2 ), even smaller fractions should be probed in 
the near future [37] . With energy loss scales 10 — 1000 Mpc for protons and 1 — 100 Mpc 
for photons both at > 10 18 eV, and source distance scales expected ~ 50 Mpc, the 
UHECR proton and photon flux collectively are perfectly suited to use as a probe of 
the local distribution of UHECR sources (discussed in the following lecture). 



LECTURE 4: MULTI-MESSENGER APPROACHES TO 

PROBLEMS 

Protons, nuclei, gamma-rays and neutrinos produced in extragalactic UHECR sources, 
such as Active Galactic Nuclei (AGN), Gamma Ray Bursts (GRBs) and Starburst Galax- 
ies, will encounter cosmic background (CB) photons during propagation through ex- 
tragalactic space. While UHE-neutrinos travel essentially without losing energy, the 
other UHE-particles may interact with the CB photons, limiting their propagation range 



to multi-Mpc scales: ~100 Mpc for protons; ~100 Mpc for nuclei; and ~10 Mpc 
for photons. This suggests the possibility of exploring cosmologically nearby sources 
(~50 Mpc) using UHECR protons, nuclei, photons and neutrinos simultaneously. 

This lecture is divided into two sections dealing with two different multi-messenger 
approaches. In the first of these sections, the proton-photon connection looks into the 
employment of the UHECR spectrum and its photon fraction to ascertain from future 
UHECR flux measurements the local distribution of UHECR sources. In this section, the 
arriving UHECR proton spectrum for different shells from a homogeneous distribution 
of sources is calculated and the photon fraction of UHECRs is obtained. The effects 
on these results due to a local source over/under-density are highlighted. In the second 
section, we focus on the UHE-neutrino flux produced directly by the UHECR source 
for the candidate sources of AGN, GRB and Starburst Galaxies. We review the main 
features of the candidate UHECR sources (GRBs, AGN, starburst galaxies). Through a 
consideration of the source opacity, dictating the relative energy losses of nucleons and 
nuclei before leaving the source region, the observation of a nuclei component in the 
arriving UHECR at Earth is used to provide strong constraints on the UHE-neutrino flux 
expected from the whole ensemble of UHECR sources. The opacity factor and neutrino 
production in these sources are discussed and the interaction rates within these sources 
obtained. We end the lecture by presenting our conclusions for the multi-messenger 
approaches discussed. 

4.1 UHECRS AND UHE-PHOTON FRACTION 

As discussed in the first lecture, UHECR protons during propagation mainly lose their 
energy through interactions with the CMB via pair creation and photo-pion production 
interactions. Thus, as discussed in the third lecture, an associated UHE-photon flux 
arises from the decay of the neutral pions, 7T° — > 77, created through the p + 7 — > p + n° 
interaction channel. These high energy photons also have their propagation range limited 
by background radiation fields (CMB and radio background). If, besides the CMB, a low 
radio background contribution is assumed, the interaction length for 10 20 eV photons is 
approximately 10 Mpc (as shown in the left-panel of Fig. 4). 

To calculate the UHECR flux at Earth, a CR source spectrum of the form 
dNcR/dEcR ~ &xp(—Ecr/ E meLX ) is adopted, with E mSLX = 10 20,5 eV, the maxi- 
mum energy CR energy attainable by the source. The value of the spectral index, a = 2, 
is motivated by non-relativistic first order Fermi acceleration (see first lecture). The 
assumed spatial distribution of CR sources is dN/dV ~ (1 +z) 3 , where z is the redshift 
of the source. Note that for the cosmologically nearby sources, any z dependence is of 
little relevance. 

To describe the contribution from regions at different distances to the UHECR flux 
(ie. to the GZK feature), five shells are taken. If L denotes a bracketed range of distances 
of the source regions from the Earth (in Mpc) and S n stands for shell n, we have: 5] for 
Le [0,3], S2 for L E [3,9], S3 for L e [9,27], S 4 for L e [27, 81], and S 5 for L e [81,243]. 

For a homogeneous distribution of sources, the arriving contribution from accelerators 
within a shell, ignoring the effects of attenuation, is proportional to the width of the 
shell. So the ratio of contributions from consecutive shells 5 n+ i to S n , is R n+ \ : R n — 



1 : 0.3. After including the effects of attenuation, it is seen that the more distant the 
source, the earlier the GZK cutoff appears. This result may be easily understood through 
an application of the analytic description for proton propagation provided in the first 
lecture. 

The contribution from the different shells (also discussed in the previous lecture) to 
the photon flux and photon fraction are shown, respectively, in the left panel of Fig. 4. For 
this result, the Clark radio background radiation and a single 7y bg interaction (ie. strong 
extragalactic magnetic fields) were assumed. From this figure it is clear that the UHE- 
photon flux is suppressed for photons coming from regions farther than tens of Mpc 
away. The UHE-photon fraction for this case, shown in the right-panel of Fig. 4, shows 
two peaks with values 10~ 3 and 10~ 2 , respectively, with the dominant contributing shell 
being 9-27 Mpc. 

To model a local over-density of UHECR accelerators, a ratio of sources from differ- 
ent shells R n +\ : R n = 1 : 0.5 and R n +i '■ R n = 1 : 0.7 was adopted. Such an over-density 
leads to a less abrupt and slightly higher in energy GZK suppression and an increase 
in the photon fraction at energies around the first peak. Similarly, to describe an under- 
density in the UHECR accelerator distribution, a local void of UHECR sources was 
introduced. Such a local void lead to a slightly lower in energy and very abrupt GZK 
suppression and a decrease in the photon fraction at energies around the first peak. For 
this case an increase in the photon fraction at energies around the second peak was also 
observed. However, such a photon component would not be detectable with present and 
(foreseeable) next generation UHECR detectors. 

The effect of diffusion on these results is to increase the path length of the protons 
relative to that of the same energy photons, reducing the flux contribution from distant 
sources and thus, in a homogeneous scenario, increasing the photon fraction component. 
In this sense the photon fraction values calculated can be considered conservative values. 



4.2 UHECR NUCLEI AND THE DIFFUSE NEUTRINO FLUX 

FROM SOURCES 



We here go through "generic" calculations for the possible neutrino flux from AGN, 
GRBs and starburst galaxies, making use of the source opacity, f K , for py interactions. 
At energies above 10 15 eV, the sources of CR remain completely unknown. However, a 
list of candidates sources may be drawn up from the Hillas criteria [38]. Acceleration 
mechanisms aside, a simple dimensional argument of the possible acceleration sites and 
the maximum obtainable energy from these sites due to the magnetic field and the size 
of the site can be considered, giving 



where E max is the maximum energy obtainable by a particle in the source, B is the 
magnetic field in the source, and 7? source is the size of the source region. Together 
with a consideration of the physics to describe the escape time of the particle and the 
energy gain rate such as that achieved through diffusive shock acceleration models under 
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Figure 6. Left: The UHECR proton flux arriving for the case of a local over-density of UHECR sources 
distribution of sources for the cases of R„+i : R„ of 0.5 and 0.7. Right: The corresponding photon 
fraction of the arriving UHECR obtained for the local over-density of UHECR sources for the cases 
of R n +\ : R„ of 0.5 and 0.7. For this plot, the label "Clark radio, single" refers to the assumptions about 
the radio background and magnetic field strengths, for which we have assumed the Clark et al. [33] 
radio background and extra-galactic magnetic fields stronger than 0. 1 nG. The injected proton spectrum 
assumed for these results had a spectral index of a=2. 



the Bohm diffusion assumption, for which the escape time ~ R 2 / (i?Larmorc) and the 
acceleration time ~ R/ (jS^c), a refined description for the maximum energy is obtained, 

£max = /^sh-Z^(5c)i? S0 urce; 

where /3 sn is the velocity, in units of the speed of light in a vacuum (c), of the shock prop- 
agating through the source. Possible accelerators, able to achieve the particle energies 
observed at the high energy end of the spectrum (~ 10 20 eV) are active Galactic nuclei 
(AGN) with R ~ 10 2 pc and B ~ G, Gamma Ray Bursts (GRB) with R ~ 10 6 pc and 
B ~ 10 3 G and starburst regions with R ~ 100 pc and B ~ 100 /lG. 

The luminosity break energy at the candidate sources is 10 44 erg s _1 for AGN, 
10 52 erg s 1 for Gamma Ray Bursts (GRBs) and 10 42 erg s 1 for starburst galax- 
ies. These numbers must then be compared with the luminosity of a UHECR accel- 
erator. With UHECR above 10 18 eV having an energy density ~10~ 8 eV cm~ 3 or 
10 54 erg Mpc 3 , the flux density of the sources, outputting for a Hubble time must 
be 10 37 erg Mpc~ 3 s _1 . With source densities ~ 10~ 5 Mpc -3 , this puts a requirement of 
L ~ 10 42 erg s _1 . Thus, all the sources considered it seems may be considered reason- 
able, physically motivated, candidates. 

AGN: In this model, UHECRs are accelerated within relativistic blobs of plasma moving 
along the AGN jet. Assuming jet Lorentz factors ~ T ~ 30, and acceleration regions of 
size /source ~ TcAt ~ 10 -2 pc, At is the typical duration of a flare, the photon density 
within the source region with photon energies at which the spectrum (EydNy/dEy) is 
maximum, is n max = 10 cm~ 3 . This leads to expected opacity factors of the source 
of f K = /source Ant. ~ 500 [39], where /source is the size of the source region and l mL ps 
1 / {nyO p y) is the interaction length of the particle. 



GRB: GRBs are brief flashes of high energy radiation that can outshine any other 
source in the sky. It is believed that the energetics associated with a GRB are caused 
by the dissipation of the kinetic energy of a relativistic expanding plasma wind called 
a "fireball". CRs may be accelerated in the fireball's shocks where the Lorentz factor 
are r ~ 300. During the post-fireball phase, due to the expansion of the spherical shock 
in the shell rest frame and the relative motion between the shell rest frame and the 
observer's frame, the size of the original object is / SO urce = T 2 cAf ~ 10 -6 pc, where 
At is the smallest observable timescale in the prompt emission. The radiation field 
spectrum is modelled according to EydNy/dEy oc EZ ^ , with /3 = 0, 1 below and above 
the break energy Eyf, r = 1 MeV , respectively [40] . The number density for photons is 
found to be n™^ ~ 10 17 cm~ 3 . This leads to expected opacity factors of the source of 
U = Wee Ant. « 0.5 [39]. 

Starburst Galaxy: These are galaxies undergoing an episode of large scale star for- 
mation. Within the starburst region of size L ~ 100 pc, OB and red supergiant stars 
emit UV radiation, whereas the dust surrounding the starburst region produces in- 
frared emission[39]. These components are describable by a blackbody and a grey- 
body respectively. In this sources, the photon number density is n r ~ 10 5 cm~ 3 , much 
less than in AGN or GRBs. This leads to expected opacity factors of the source of 
h = WceAnt. « 10- 4 [39]. 

Protons interact with the source ambient radiation to photo-produce pions when the 
photon energy in the proton's rest frame is above Ey-th ~ 145 MeV. The decay of the 
secondary neutrons and charged pions produce a flux of neutrinos (as discussed in lecture 
2). 

Both for AGN and GRB, charged pion decay is the dominant contribution to the 
cosmic neutrino flux, which at its peak value is around 75% of the Waxman-Bahcall 
bound. For the Starburst model, due to the small values of f n , the neutrino flux is 4 
orders of magnitude below the Waxman-Bahcall limit. 

If UHECR sources accelerate heavy nuclei instead of protons, the resulting high 
energy neutrino spectrum will be suppressed compared to an all-proton composition 
at acceleration. At the source, a nucleus interacts with the radiation fields photo- 
disintegrating into its constituent nucleons, which may then themselves produce neu- 
trinos through photo-pion interactions. Since the density of the radiation fields at the 
source determines the intensity of the associated neutrino flux, if most of the nuclei are 
broken into their nucleons, the flux of neutrinos will be very similar to that predicted for 
accelerated protons. On the other hand, if the accelerated nuclei leave the source without 
interacting, the cosmic neutrino flux will be dramatically suppressed. 

We can apply the proton interaction rate (as discussed in the first lecture) 

M pr +A pr )/2T P dn Y (£ Y ) 

Rpy^Opy d£y^^, (19) 



£ P7 -A py )/2T p d£ y 



to photo-disintegration reactions 



A£, Fey +h Fe7 )/2Tp driy(£y) 

Rfey^Ofey / d£y ; . (20) 

J{^ Fey -A Pey )/2T p d£y 



where driy/d£y is the number density spectrum of ambient photons, o p y ~ 0.5 mb, 
% pr « 310 MeV, A pr « 100 MeV, o FeJ « 80 mb, £ Fer « 18 MeV and A Fey « 5 MeV. 
This gives a ratio between photo pion and photo-disintegration rates of 

W r p) « — ^r( 15 r p) = I60^ y (i5 r p ) . (21) 

Complete photo-disintegration only occurs in GRBs for Ep e > 10 17 eV, and in AGN, for 
Ef e > 10 20 eV. The photo-disintegration at Starburst regions is negligible. As expected 
from this expression, in all cases the photo-disintegration curve followed the shape of 
the proton interaction rate, with the shape of both curves being dictated by the shape of 
the radiation field spectrum. 

4.3 CONCLUSIONS 

In the first part of this lecture, the detection of the photon fraction component of UHECR 
spectrum has been demonstrated to be a powerful diagnostic tool for both verifying the 
origin of the suppression feature observed of the in the UHECR spectrum as well as 
being usable, in conjunction with the UHECR cut-off feature observed, in determining 
the local distribution of UHECR sources. In the second part of the lecture, generic 
calculation for the diffuse UHE-neutrino spectra from candidate sources of UHECR 
(protons) have been gone through. Through the assumption, in both GRB and AGN 
models, for near unity neutrino production opacities (/^), complete disintegration of 
UHECR nuclei in these objects has been demonstrated to be expected fpr energies 
above 10 17 eV and 10 20 eV respectively. A consistency check of these opacities factors, 
whose values dictate the probability of UHECR protons interacting and producing UHE- 
neutrinos in the sources before escaping, with the presence of Fe in the arriving UHECRs 
demands that much smaller neutrino fluxes from these sources should be expected than 
obtained from present calculations which ignore the information carried by the existence 
of UHECR nuclei. 
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